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Abstract 

Energy harvesting is increasingly gaining importance as a means to charge battery powered devices 
such as sensor nodes. Efficient transmission strategies must be developed for Wireless Energy Harvesting 
Nodes (WEHNs) that take into account both the availability of energy and data in the node. We consider 
a scenario where data and energy packets arrive to the node where the time instants and amounts of the 
packets are known (offline approach). In this paper, the best data transmission strategy is found for a 
finite battery capacity WEHN that has to fulfill some Quality of Service (QoS) constraints, as well as the 
energy and data causality constraints. As a result of our analysis, we can state that losing energy due to 
overflows of the battery is inefficient unless there is no more data to transmit and that the problem may 
not have a feasible solution. Finally, an algorithm that computes the data transmission curve minimizing 
the total transmission time that satisfies the aforementioned constraints has been developed. 
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I. Introduction 

Battery autonomy is limiting the operational lifetime of battery powered wireless nodes due 
to the increasing computational capabilities, which follow from the reduction of transistor's size, 
and the consequent increment of the node throughput requirements. As battery capacity is not 
growing as fast as the energetic demand of battery powered wireless nodes, e.g., handheld devices 
or sensor nodes, it is essential to make an efficient use of the available energy. 

The energy consumption to transmit a certain amount of data can be reduced by increasing the 
transmission time or, equivalently, reducing the transmitted power or rate at which the information 
is sent [1]. However, there exists a trade-off between energy and transmission delay as in many 
real scenarios, the data flow must guarantee certain levels of performance. For instance, the 
network may require that data is transmitted with a certain minimum rate, quality level or with 
a maximum delay. These requirements, usually referred to as QoS constraints, become crucial 
in real-time applications that are sensitive to the network delay. 

The trade-off between energy and transmission delay was studied in different works [1]- 
[5], where the objective was to minimize energy consumption while satisfying QoS constraints. 
Specially remarkable is the work by Zafer et al. in [1], as they proposed a novel solution that 
used the concept of cumulative curves, which allows a graphical and simple visualization of the 
solution. 

Within this context, energy harvesting, i.e., the process of collecting energy from the envi- 
ronment by different means, e.g., a solar cell or a piezoelectric generator [6], [7], is becoming 
an appealing solution to enlarge the autonomy of battery powered devices. Energy harvesting 
has opened a new research paradigm on the design of energy-efficient transmission strategies as 
the transmission strategies for non-harvesting nodes are no longer optimal. Many works, e.g., 
[8]-[16], model the energy harvesting process as a set of energy packets arriving to the node 
at different time instants and with different amounts of energy. There exist two well established 
approaches for the design of optimal transmission strategies, namely, online and offline. The 
online approach assumes that the node only has some statistical knowledge of the dynamics of 
the energy harvesting process. The offline approach assumes that the node has full knowledge of 
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the amount and arrival time of each energy packet, which is an idealistic situation that provides 
analytical and intuitive solutions that can be used to gain insight for the later design of the online 
transmission strategy. 

By using this packetized model for the energy harvesting process, [8] considered dynamic 
data packet arrivals and found the transmission strategy that minimizes the delivery time of all 
data packets, however, the authors of [8] assumed an infinite battery capacity. In [9], a node with 
finite battery capacity was studied, however, considering that all the data packets are available 
from the beginning of the transmission. This assumption significantly simplifies the setup as 
losing energy due to battery overflows is clearly suboptimal because there is always data to 
be transmitted. References [10] and [11] used convex optimization and dynamic programming 
to find the best rate scheduling or power allocation strategy under the presence of fading. The 
minimization of the transmission completion time for a WEHN operating in a broadcast link 
was considered in [12]-[14]. References [12] and [13] assumed infinite battery capacity, and, 
the authors of [14] found the rate scheduling policy of the finite battery capacity case. However, 
these works assumed that the data to be transmitted to each of the users is available from the 
beginning at the transmitter, similarly as in the point-to-point case in [9]. In [15], the coding 
problem was studied from an information theory perspective for an energy harvesting node. 
Finally, the stability and delay of the data queue were considered in [16] to derive the optimal 
transmission policies. 

To the best of our knowledge, there is no work that considers a WEHN with finite battery 
capacity and dynamical data and energy arrivals. In a prior work [17], we studied this problem 
and showed that as far as battery does not overflow, constant rate transmission is the strategy 
that requires less energy to transmit a certain amount of data. However, if battery overflows, 
transmitting at constant rate is no longer optimal, but the optimal strategy increases the rate 
before the overflow until either there is no overflow or until all the data is transmitted. In this 
paper, we extend the work in [17] to take into account new QoS requirements of the different 
data packets, which was not considered in [17]. QoS constraints substantially complicate the 
problem as we have constraints both in the data and energy domains. In this paper, we propose 
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a framework to map the constraints of the energy domain to the data domain that allows us to 
adapt the calculus approach proposed in [1], which did not take into account energy harvesting 
at the transmitter, to the energy harvesting scenario. Therefore, the main contributions of this 
paper are: 

• Studying the impact of the QoS constraint in the transmission strategy that minimizes the 
transmission completion time, which to the best of our knowledge has not been previously 
considered in the literature. 

• Showing that, due to the QoS constraints, there may be situations in which no feasible 
transmission strategy exists. Such situations are analytically characterized in Lemma 4. 

• Showing that if the optimal transmission strategy exists, the optimal cumulative data depar- 
ture curve is a peace-wise linear function where battery overflows are only produced when 
the data buffer is empty. 

• Developing an algorithm that either computes the optimal transmission strategy or concludes 
that there is no feasible solution and analytically showing its optimality. 

In Section II, the problem is mathematically formulated by using the concept of cumulative 
curves, which is introduced in [1], that allows an appealing visualization of the solution. The 
solution is characterized in Section III. Section IV, presents the developed iterative algorithm 
that is able to compute the optimal solution by using the characterization performed in Section 
III. Finally, the conclusions of the work are presented in Section VI. 

II. Problem formulation 

We consider a node with a finite battery capacity, Cmax, that has to transmit N data packets 
by using at the most the J energy packets that it harvests over time while satisfying some QoS 
requirements (see Figure 1). We want to find the power allocation/rate scheduling strategy' that 
minimizes the transmission completion time, T. 

We assume that the time instants at which the data and energy packets arrive to the node 
and their size (bits or Joules) are known from beforehand {offline approach). Hence, it is known 

'observe that fixing the transmission power or the rate is equivalent as it will be shown next. 
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Figure 1. System Model. 




Figure 2. Summary of the different events considered in the time domain, namely, energy arrivals, data arrivals, and quality 
of service requirements. The figure represents the deadline QoS constraint, where 9 is the maximum allowed delay for all the 
packets. 

that at the time instant di > seconds the i-th data packet arrives containing Dj bits, with 
i = ... N — 1. Similarly, the j-th energy packet arrives at the instant ej > seconds and a 
total of Ej Joules are harvested, with j = ... J — 1 (see Figure 2). Without loss of generality, 
the first energy arrival is produced at eo = and contains the initial battery of the node Eq. 
To describe our model we present the following definitions that are summarized in Table I: 
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Definition 1 (Data Departure Curve). A data departure curve D{t) > 0, t > 0, is the total 
number of bits that have been cumulatively transmitted by the node in the time interval [0,t]. 

Definition 2 (Energy Expenditure Curve). An energy expenditure curve E{t) > 0, t > 0, is the 
energy in Joules that has been cumulatively consumed by the node in the time interval [0, t]. 

Let us consider a static or slow-fading channel where the power-rate function g{-), i.e., the 
function that, at any given time instant t, relates the transmitted power, P(t), with the rate, 
r(t), according to P{t) = g{r{t)). As in [1] and [8], we make the common assumption that 
the function g{-) is time-invariant, convex, monotonically increasing, and only depends on r{t). 
Note that the instantaneous rate, r(t), can be expressed as the derivative with respect to t of the 
data departure curve, i.e., r{t) = Similarly, the transmitted power is P{t) = Then, 

the energy expenditure curve can be obtained from the data departure curve as follows^: 



Observe that the magnitudes D(t), E(t), r(t), and P{t) are unambiguously related by (1) and 
g{-). Therefore, given the initial states E{0) = and -D(O) = 0, the design of the system to be 
optimized can be described by any of these magnitudes. 

Definition 3 (Battery). The battery of the node Bit) is the amount of energy that the node has 
available at a given time instant t. We consider a battery with finite capacity Cmax- Thus, B{t) 
must satisfy that < B{t) < Cmax, Vt > 0. 

Due to the limited battery capacity, at the j-th energy arrival, some part of the harvested 
energy Ej may be lost. This lost energy is denoted as the j-th battery overflow, i.e., Oj = 
{Ej — Cmax + B{ej)}^, with {x}^ = max{0, x}. Observe that battery overflows depend on the 
chosen energy expenditure curve E{t) and guarantee that the battery level will never be above 
the battery capacity, B{e^) = B{eJ) + Ej — Oj < Cmax- 

^Without loss of generality, we can assume that r{t) is right-continuous. 




(1) 
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As battery overflows depend on the chosen E{t), their value Oj cannot be computed until 
the energy expenditure curve E{t) is fixed Vt < ej. Loosely speaking, one can expect that 
the optimal solution uses efficiently the harvested energy to transmit the available data and, at 
the same time, tends to minimize the total overflow of the battery and, thus, maximizes the 
accumulated energy stored in the battery. In other words, if overflows are minimized, the node 
will be able to use more energy, as, otherwise, the energy of the overflows is lost. 

In the following lines, we define the accumulated battery, a concept introduced in this work 
that allows us to characterize the optimal solution when having finite battery capacity constraint. 
Let denote the last time instant up to which D{t) (or, equivalently, E{t)) is known. 

Definition 4 (Accumulated Battery). The accumulated battery BA{t]tx) is a real measure of 
the accumulated energy stored in the battery for t G [0, t.j] and it is the maximum possible 
accumulated energy in the battery for t G {t,j.,oo) (assuming that no overflows are produced 
for t G (tx, oo)), i.e., BA{t] t^) = Y.j=o {-OjH{t^ - ej) + Ej) H(t - ej), where H{-) is the 
Heaviside function. 

Observe that, for t G [0,^^^^], BA{t]tx), represents the real measure of the energy accumulated 
in the battery during (0, t), because, for these time instants, overflows are known and taken into 
account. Alternatively, for t G {tx.oo), battery overflows are unknown after tx and BA{t;tx) 
models the best case scenario where the node is able to store in the battery all the harvested 
energy in the interval (tx, oo). Moreover, observe that BAit; tx) = BA{t; ty) for any tx, ty > t. 

At every time instant, the energy stored in the battery is B(t) = BA{t]t) — E(t). Note that 
BA{t;t) takes into account the actual net incoming energy in the battery, whereas, E{t) is the 
net outgoing energy. Thus, their difference results in the energy stored in the battery. 

Definition 5 (Minimum Energy Expenditure). The minimum energy expenditure, Emin{t]tx), 

is the smallest amount of energy that the node must have cumulatively spent at time t > tx 
such that no overflow of the battery is produced in the interval (tx,t], i.e., Emin{t]tx) = 
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Definition 6 (Accumulated Data). The accumulated data Dy^{t) is the sum of data that has 
arrived at the node during the time interval [0,t), i.e., DA{t) = Xlilo^ DiH{t — di). 

Different QoS constraints can be considered by mapping the constraint into an appropriate 
minimum data departure curve, which was introduced in [1], and is defined as follows: 

Definition 7 (Minimum Data Departure). The minimum data departure, Dqosit), is the smallest 
amount of data that the node must have cumulatively transmitted at time t such that the QoS 
constraint is satisfied. 

The rest of the paper considers a general Dqosit) that is a non-negative staircase function 
where changes are produced at time instants qk vvith increments of Qk bits for k = . . . K — 1, 
i.e., DQos(t) = Yl!k=Q QkH{t — Qk). From now on, the instants qk are called quality requirement 
events. Thus, three kind of events are considered, namely, data arrival, energy arrival, and quality 
requirement events, as summarized in Figure 2. Let us briefly describe two examples, introduced 
in [1], of how to map a certain QoS constraint to the minimum data departure curve: 

• Deadline Constraint: This constraint considers that the maximum permissible delay for the 
transmission of a certain data packet, D^, is 6^ seconds. Then, Dgosit) is a piecewise 
constant function that changes at instants Qk = dk + 9k with an increment of Dk. As a 
specific case, we can consider that the allowed delay for all the packets is the same, i.e., 
9k = 9, \fk. Then, the minimum data departure is given by DQos{t) = DA{t — 9). 

• Buffer constraint: A limited data queue of size (3 is considered. Hence, in order not to loose 
any incoming data, the minimum data departure must be Dgosit) = {-^^(t) — Z?}"*". 

Observe that depending on the chosen QoS, it is likely that the instants qk are equal to for 
some values of k and i, e.g., for the buffer-type constraint. Hence, we consider that two different 
types of events can be produced simultaneously at the same time instant.^ 

Our goal is to find the data departure curve, D(t), that minimizes the transmission completion 
time T of the N data packets, i.e., D(T) = J2f=o^ while satisfying the following conditions: 

^Because of this possible simultaneity between events, the problem considered in this paper is much more difficult than the 
one considered in [17]. 
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(i.) Energy Causality Constraint (ECC): energy must be harvested before it is used by the node 
or, which is the same, the battery level in the node must be greater or equal to zero, (ii.) Data 
Causality Constraint (DCC): it is not possible to transmit more bits than the ones that have 
arrived to the node. (Hi.) QoS constraint: at time t, a minimum amount of data DQos{t) has 
to be transmitted in order to preserve the link quality of service. Moreover, given two data 
departure curves with the same completion time, the one that requires less energy is always 
preferred. From all that has been said above, the problem can be expressed as follows: 

min T (2) 

D{t) 

s.t. E{t) < BA{t]t), 

DQos{t) < D{t) < DA{t), 

N-l 

We want to remark the two main difficulties of the problem presented in (2). First, the integral 
relation among the data and energy domains through (1). Second, the fact that neither T nor 
BA{t;t) are known from beforehand, due to their dependence on D{t). Consequently, both T 
and t) will be found along with the solution to the problem. 

This problem is graphically represented in Figure 3, where the figures at the top and bottom 
stand for the energy and data domains, respectively. The ECC is represented by the solid line 
at the top figure, whereas DCC and QoS constraints are depicted by the dot-dashed and dashed 
lines in the figure at the bottom, respectively. The dotted line in the energy domain represents 
the minimum energy expenditure curve, however, as it can be seen from (2), it is not a constraint 
of the problem. Hence, D(t) and its associated E{t) must lie within the blank region of the data 
and energy domains, respectively, in order to be a valid solution. Three different data departure 
curves (A, B, and C) and their associated energy expenditure curves are shown. The curve A is 
not valid since it breaks the ECC. The curve B is valid in spite of having an overflow of the 
battery at 63, which is, in general, a suboptimal strategy as we will show later. Finally, the curve 
C is not valid because it does not satisfy the QoS constraint. 
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Figure 3. Illustrative representation of the problem presented in (2). 

We want to remark that Figure 3 is just an illustrative representation of the problem. As 
mentioned before, both the accumulated battery and the minimum energy expenditure depend 
on the selected transmission strategy. Note that the figure shows BA{t; 0) as ECC instead of the 
real ECC BAit; t). In some manner, with BA{t] 0), we are showing the battery accumulated that 
is obtained when D{t) does not produce any battery overflow. In case that overflows cannot be 
avoided, it would be necessary to subtract the amount of the overflow to the shown accumulated 
battery and minimum energy expenditure from that time instant onward. Thus, we can only show 
the real graphical representation of the problem once we have fixed the solution. 

Note that (2) is not a convex optimization problem and that its conversion into a convex 
problem is not straightforward. Thus, we cannot directly solve (2). Alternatively, in next section, 
we model the properties that the optimal solution must satisfy, which will allow us to construct 
the optimal data departure curve. 
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III. Properties of the optimal solution 

As pointed out previously, in this section we will characterize the optimal data departure 
curve, D*{t), and its associated energy expenditure curve, E*(t), for the problem (2): 

Problem 1 (Transmission Without Events). Consider that the optimal departure curve is known 
up to ti and that we want to characterize its optimal behaviour in the time interval (^1,^2) 
where there are no changes in Dj^{t), and Dgosit) ■ We also consider that the data 

departure curve at the boundary of the intervals is D(ti) and D{t2), respectively, and that these 
two points satisfy the data, energy, and QoS constraints. 

Lemma 1. In Problem 1, D*{t) is a straight line where the slope, or, equivalently, the trans- 
mission rate, is constant and equal to rit) = "^^^^^if ^^^^ Vt G (^1,^2)- 

Proof: The proof follows from the integral version of Jensen's inequality in a similar way 
to the BT-problem in [1]. ■ 

Corollary 1. Lemma 1 implies that D*{t) is a piece-wise linear function such that its slope, 
which is equivalent to the transmission rate, can only change either at ej, di or q^. 

From the previous lemma and corollary, it follows that constant rate transmission saves energy 
due to the convexity of the power-rate function g{-). However, constant rate transmission is not 
optimal when a battery overflow is produced because the energy saved due to constant rate 
transmission is lower than the energy lost in the overflow. Consequently, the optimal solution 
increases the rate before the overflow until either there is no overflow or the data buffer is empty, 
as shown in the following lemma and its subsequent proof: 

Lemma 2. Under the optimal policy, battery overflows may only be produced when there is no 
data to be transmitted. 

Proof: A similar proof was given in [17] for the QoS relaxed problem. Note that Lemma 
2 states that if an overflow is produced, the data departure curve must reach its upper bound, 
i.e., D(t) = DA{t). Thus, the insertion of the QoS constraint, which is a lower bound on D(t), 
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does not affect the validity of the proof given in [17]. Therefore, Lemma 2 can be proved by 
following the approach in [17, Appendix A]. ■ 
The following lemma states that by the end of the transmission the battery must be empty, 
otherwise, transmission could have been finished earlier by transmitting at a higher rate. 

Lemma 3. The optimal solution must satisfy that, at the instant T at which all the data has been 
transmitted, the energy expenditure is equal to the accumulated battery, i.e., E*{T) = Ba(T; T). 

Proof: Similar to the proof of Lemma 5 in [9] . ■ 
In the remainder of the paper, the term epoch denotes each of the time intervals of D*(t) at 
which transmission is done at constant rate. We define M as the total number of epochs of the 
optimal solution, i.e., the number of linear pieces of D*{t). Note that M is unknown a priori. 
Consequently, solving the problem in (2) is equivalent to determining the rate and length irn 
of each epoch, i.e., {rm,^m}m=o^^- To do so, we have developed an iterative algorithm that, at 
the m-th iteration, determines the rate and duration of the m-th epoch, i.e., {r^, ^m}- We denote 
Tm as the instant at which the m-th epoch begins, i.e., = YH^^Iq f-p- The algorithm ends when 
all data has been efficiently transmitted. With this, M, T, BA{t;T), and D*(t) are found. 

To simplify the complexity of our algorithm, at the beginning of the m-th iteration"^, the origin 
of coordinates is moved to the point (rm , D*{Tm))- To be coherent with the vertical and horizontal 
displacement of the origin of coordinates, the data and energy constraints in (2) must be vertically 
rescaled by D*{Tm) or E*{Tm), respectively, and temporally displaced by Tm- In the remainder of 
the paper, a super-index (m) above a variable (B^J^\t), D^^\t), DQl^g{t), and E^J^^{t)) denotes 
that it is the rescaled version at the m-th iteration^, e.g., = BA{t + r^; Tm) — E*{Tm)^. 

From the structure of the problem in (2), it is easy to expect that it may not have a feasible 
solution whenever the node has to fulfill very tight QoS requirements, while, at the same time, 
it does not harvest enough energy to transmit all the required data. The following lemma is 

■^Note that D*{t) is known in [0, Tm]- 

^The relations among the iteration specific and general versions of the variables are given in Table I. 

*We have dropped the second argument in B^™' (t) and -B^^ (t) since within the m-th iteration, the second argument, which 
denotes the last instant at which the solution is known, is always r„i. 
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checked at every iteration of our proposed algorithm to determine whether the problem in (2) 
has a feasible solution or not. 

Lemma 4. The problem (2) does not have a feasible solution whenever 

^QosK(lk )>(lk 9 [^A {Qk )/% ) > (^) 

for some quality requirement event q^^^ G (0,T(™))^. 

Proof: Whenever we encounter that the problem does not have a feasible solution is because, 
at some quality requirement event it is not possible to fulfill all the constraints. Let D{t) be 
the data departure curve that transmits the maximum amount of data in the interval (0, g^,'"^), i.e., 
D(t) = {^B'^A\'lk^'') / 'Ik^^^ Ik"''- Note that this curve has constant rate, empties the battery 
at ql"^\ and that the constrains are not necessarily satisfied. If, at some g^™^ G (0,T(™)), the 
QoS constraint requires more than g^^ (^-^A^^l^'i"^^)/^!™^) Q'i™^ ^i^^ transmitted, then the 

problem does not have a feasible solution. ■ 
If there is no G {0,T^"^'>) that satisfies Lemma 4, then the problem still may have a 
feasible solution and at least another epoch, i.e., {rm,(^m}, can be determined. In this context, 
in next subsection, we model how the rate changes must be produced in order to be optimal. 

A. Constraints mapping into the data domain for a given epoch 

Within an algorithm iteration, the ECC can be mapped to the data domain, hence, allowing us 
to merge both constraints to the most restrictive constraint. Let us consider that the algorithm is 
at the beginning of the m-th iteration, i.e., the optimal solution is known up to r^, where the rate 
Tm and length im of the m-th epoch must be determined. Given that transmission must be done at 
constant rate/power, the maximum amount of data that can be transmitted at a certain time instant 
ty due to the energy causality constraint is g'^{py)ty, where Py = B^^\ty)/ty.^ With this, at the 

^Where x'"*^ = X — Tm is rescaled version of some temporal variable X at the m-th iteration (see Table I). 

*Note that for the constraints mapping it is not necessary that Py satisfies the constraints. The constraints fulfillment is enforced 
by the algorithm that computes the optimal solution, which is explained in Section IV. 
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instant ty, the upper-bound on the energy expenditure curve has been mapped to an upper-bound 
on the data departure curve, as shown in Figure 4 for the instants ti and t2- By applying this 
procedure at all time instants ty G {0,T^"^^), we can map the whole upper-bound in the energy 
domain to an upper-bound in the data domain, which we denote by D^^^{t) = g^^{B^^\t)/t)t 
and call actual mapping. However, doing this computation for each time instant has a high 
computational cost. In summary, if a data departure curve that transmits at constant rate, i.e., 
= r„t, satisfies that D^'^^t) < D^^}{t),\/t, then it also satisfies the ECCs. 

The cost associated with the computation of D^^^{t) can by reduced by noting that it is 
suboptimal that D^^^^t) reaches the actual mapping at any time instant ti that is not an event, 
i.e., ti ^ d["^\ ti 7^ and ti ^ (lk"\ Vz,j, A; . This is clearly seen in Figure 4. Observe 
that if = D^^}{ti), then the battery is empty at ti, i.e., E^'^^ti) = B^^\ti), which 

follows from the definition of the actual mapping. Consequently, in order to satisfy ECCs, the 
rate at t]^ must be zero as no energy arrival is produced at ti. From Corollary 1, we know that 
this rate change is suboptimal and, therefore, the data departure curve can only reach the actual 
mapping in some event. Thus, to reduce the computational complexity of the actual mapping, 
we can compute the value of the mapping only at the aforementioned events and assign a 
constant value in the interval between events. In the rest of the paper, we refer to this mapping 
as ejfective mapping, i.e., D^^^(t). Note that the ejfective mapping is an upper bound of the 
actual mapping. However, we want to remark that by using the effective mapping, we are not 
relaxing the constraints of the problem since both mappings are equal at the time instants where 
the optimal solution for the data departure curve coincides with the actual mapping. 

A similar approach can be done to map the minimum energy expenditure curve El^l{t) to 
the data domain. In this case, since overflows may only be produced at energy arrival events, it 
is only necessary to map the lower bound in the expended energy to the data domain at these 
time instants. For the rest of time instants, the time intervals between energy arrivals, a constant 
value is assigned without loss of generality, hence, obtaining D^]^^^^{t). 

Figure 5 shows a representation of the problem once the energy constraint is mapped to the 
data domain. Now the problem is simplified, since data and energy constraints can be merged in 
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(m) 



P2 = B!r'(t2)/t2 



t2 




Figure 4. Mapping of thie energy constraint to the data domain. The plot in the bottom shows the suboptimality to reach 
Dg^ (t) at any time instant that is not an event. This implies that in practice the effective mapping D^j^ (t) can be used as the 
mapping of the energy constraint to the data domain. 



a single constraint that, at every time instant, is the most restrictive of the two constraints, i.e.. 



^LlW=min{Di™)(t),<)(t)}. 



(4) 



Similarly, the lower constraint is 



(m) 



Jm) 



(5) 



Note that D^Jx{t) and D^l^{t) are only valid within an algorithm iteration and that it can 
occur that Dlnax{t) < D^l^{t). This happens in the following situations: (i.) The node has to 
transmit a certain amount of data in order not to overflow the battery, however, this data is not 
still available, i.e, D^^_^{t) > D^^\t). (ii.) The node has to transmit a certain amount of data in 
order to satisfy the QoS constraint, however, it does not have enough energy available to do so, 
i.e., D^^^{t) < DQ^g{t). Note that the situation (ii.) occurs when the problem does not have a 
feasible solution. As mentioned before, the aim of this section is to model how rate changes are 




Figure 5. Mapping of the ECC and minimum energy expenditure to the data domain. 



produced when the problem indeed has a solution (at least up to the current algorithm iteration) 
and, hence, we will focus on situation (i.) where an overflow of the battery is produced. 

Let us define the sets of time instants at which Dmax{t) and D^l^{t) have discontinuities 
as ztL = {t I DlTlit") ^ DtUt+)} and Z^^l = {t \ D^^lit^) ^ D^Zlit^)}, respectively. 
Remember that, due to Corollary 1, we know that D*{t) is constant between events defined 
according to and Z^^. In Lemmas 5, 6, and 7, we model the behavior of the optimal 

solution when the rate changes at a time instant where a single event is produced. Similarly, 
Lemmas 8 and 9 describe the behavior of the optimal solution when the rate changes at a time 
instant where two events are produced. The proofs of these lemmas are given in Appendix A. 

Lemma 5. If a rate change is produced at a certain time instant £m such that im G Z^x and 
£m ^ ^l^il, then D*^"'\£rn) = Dlnax{im) ^^^^ ^^^^ incrcascs, Vm < Tm+l. 



Lemma 6. If a rate change is produced at the time instant im such that im ^ Zinax, im ^ 
and -D^^l(£m) > D^^l{i:^), then D*'''"\tm) = ^'^^ the rate decreases, r„ > Vm+i- 



Lemma 1. If a rate change is produced at a certain time instant im. such that im ^ zim^x, 
im e and DlnJx{im) < -D^^(^m), then an overflow of the battery is produced at im, 
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D*''"'\irn) = D^xiim) ci^d the rate is zero until the next data arrival event, r^+i = 0. 

Lemma 8. If a rate change is produced at a certain time instant irn such that im G 

and im e and Dl^i{i-) > D^ZliO' then either = ^^al(C) and the rate 

increases, or D*''"^\im) = D^'lli^) and the rate decreases. 

Lemma 9. If a rate change is produced at a certain time instant im such that im G and 
im e Z^l, and DlnJx{im) < D^l^{im), then an overflow of the battery is produced at im, 
D*^"^\tm) = Dlni^x{im) ^^e rate can either increase or decrease. 

By using these lemmas, we are able to construct an algorithm, which is presented in next 
section, that iteratively finds the optimal solution, or concludes that there is no feasible solution. 

IV. Optimal data departure curve construction 

In this section, we describe the developed algorithm that is able to either construct D*(t) 
or, alternatively, conclude that the problem in (2) does not have a feasible solution. As stated 
in Corollary 1, the optimal data departure curve is a piece- wise linear function. As previously 
explained, the developed algorithm follows an iterative process where, at the m-th iteration, the 
duration, im, and rate, r^, of an epoch are determined. We will focus on the explanation of the 
m-th iteration since all the other iterations follow the same approach. 

As shown in Figure 6, the algorithm is composed by three main blocks. The first block, 
named check Solution, determines the existence of solution in the current iteration by checking 
the condition in Lemma 4. If the problem does not have a solution, the algorithm ends. Otherwise, 
the algorithm proceeds to the subsequent blocks to determine the rate and length of the epoch. 

The second block, named checkFinish, checks whether it is possible to transmit all the 
remaining data by using all the available energy at constant rate. This block is necessary because 
T^^'> is unknown. In case it is possible to finish in a single epoch, it is obtained that T*^™^ = im 
and M = m + 1. Then, the solution to problem (2) has been found and can be computed as 

M-l fl_T-\ 
D\t) = rm{t - r™) n ( — ^ ) + rmimH{t -Tm- im), (6) 
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m 



checkSolution 




yes 

m + + I 



checkFinish — Yes- 



M 



= m + 1 I 



No 



getEpoch 



Mode 
{minT or min Energy) 



Figure 6. Block diagram of the iterative algoritlim. 



where n(-) is the unit pulse in the interval [0, 1]. The domain of D*(t) is [0,T = J2m=o ^m]- 

Otherwise, if it is not possible to finish in a single epoch, the block checkFinish returns the 
mode (minT or minEnergy) to be used by the third block, which we name getEpoch, to find 
the epoch rate and length (r^ and irn) that fulfill Lemmas 5-9. The mode minT is used when 
the node already has enough energy to finish transmission, whereas, the mode minEnergy is 
used when the node is still not able to finish transmission at any rate and, hence, the objective 
is to save as much energy as possible for the end of the transmission. An extended explanation 
of the inner behavior of each of these blocks is given in Appendix B. 

Once r„ and £m are determined, the origin of coordinates is moved to the point {£m, rmim) and 
the variables are prepared for the new iteration. In the data domain, the iteration transmitted bits 
D'^*'"^^) = r^i^ are subtracted from D^^\t) and ^'^^^(t), for instance, D^^^^\t) = D^^\t+ 
im) - D*^"'\im). Similarly, in the energy domain, the expended energy E*'^"'\im) = g{rm)im 
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is subtracted from B^^\t) and E^l^{t). Moreover, in case that transmitting at produces 
a battery overflow at time instant the amount of energy lost due to the overflow is also 
subtracted from these variables. Finally, the mapping to the data domain is recalculated for the 
iteration m + 1 and the whole procedure starts again to determine r^+i and Im+i- 

There are two possible reasons for which the algorithm ends: (i.) At some iteration. Lemma 4 
is satisfied and, hence, the problem does not have solution, (ii) All the data has been transmitted 
and the optimal data departure curve has been obtained as given in (6). 

The algorithm optimality is summarized in the following theorem and its subsequent proof: 

Theorem 1. The algorithm presented in this section constructs the optimal data departure curve, 
D*(t), for the problem (2). 

Proof: See appendix B-D. ■ 

V. Results 

To the best of our knowledge there is no other algorithm in the literature that considers 
altogether the ECC, DCC, the QoS constraint and the finite battery capacity. Therefore, in order 
to get some insights on the gain obtained with our proposed solution, we have developed a 
suboptimal ad-hoc strategy, namely, the Empty Buffers Strategy (EBS), that tries to empty the 
buffers as soon as possible, i.e., it looks for the time instant at which the next arrival (energy 
or data) is produced and tries to transmit at a constant rate so that the corresponding buffer is 
emptied by the time of the corresponding arrival. 

The left y-axis of Figure 7 compares the normalized mean minimum T along a total of 1000 
iterations, where at each iteration the data and energy arrivals are randomly generated following 
a uniform distribution. The amount of energy in each of the packets is normalized according to 
the total harvested energy which varies along the a;-axis. The right y-axis shows the percentage 
of iterations in which there exists a feasible solution to (2). As shown in Figure 7, our proposed 
optimal algorithm substantially reduces the mean minimum T. If feasible solutions exist the 
optimal algorithm finds the one that minimizes T, however, for some arrival profiles, the EBS 
is not able to find any feasible solution in spite of its existence. 
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Figure 7. Performance comparison of our proposed optimal algorithm with respect to the suboptimal BBS in terms of mean 
minimum completion time and percentage of solutions. The lines marked with rectangles and diamonds refer to the left axis, 
whereas the lines marked with circles and asterisks refer to the rights axis. 



Note that the EBS performs worse than the optimal strategy because it changes the rate at 
every packet arrival without checking whether constant power transmission is feasible between 
two arrivals, which would consume less energy. As a result of this extra energy consumption, 
the mean minimum T is higher and the probability of finding a feasible solution decreases. 

VI. Conclusion 

The presence of energy harvesters in wireless nodes implies a loss of optimality of the 
traditional transmission policies such as the well-known water- filling strategy [18]. In this paper, 
the optimal transmission strategy has been obtained for wireless energy harvesting nodes with 
finite battery capacity that, additionally, have to fulfill some QoS constraint. Hence, we have 
contributed into the decrease of the required transmission completion time and, thus, increased 
the overall efficiency in the use of the harvested energy. Moreover, in more technical terms, we 
have seen that, as far as battery does not overflow, constant rate transmission is the strategy that 
requires less energy to transmit a certain amount of data. However, if indeed battery overflows, 
transmitting at constant rate is not optimal anymore, but the optimal strategy increases the rate 
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before the overflow until either there is no overflow, or the DCC is reached (i.e., there is no 
more data to transmit). We have seen that the existence of the optimal solution depends both 
on the dynamics of the harvesting process and on the required QoS. According to this, we have 
developed an algorithm that is able to determine whether the problem has a solution or not and, 
in case of having a solution, determines the optimal data transmission strategy. 

Appendix A 

Proof of Lemma 5 

The proof of the lemma is divided in two parts. We first show that if a rate change is produced 
at Ira, then D*''"^\£rn) = Dlnax{im) ^i^^, aftcrwards, we show that the rate must increase. 

Part 1. Let us assume that a rate change is produced at such that D'^'^\£rn) < Dinax{(m)- 
We will show by contradiction that this cannot be optimal. Let us consider the time interval 
{irn — ^Am + e) with 6 being positive. Note that if we select a sufficiently small e, we can 

find a straight line with rate r = ^ '^^ — — -, which still satisfies the constraints, that 

transmits the same amount of data while having less energy expenditure. Hence, we have proved 
that if the rate changes at i^, then D*^'^\i^) = D^li^^r.)- 

Part 2. Now we prove that when the rate changes at it must increase. The procedure is 
the same as the one in the first part of the proof. We start by assuming that a rate decrease is 
optimal and then, we see that it leads to a logical contradiction. We denote and rm+i the 
rates before and after respectively, where > Vm+i- We consider the same time interval. 

Ty(m) ff -\-f,'] nC"^) (f 

In this case, we can also find a straight line whose slope is r = ^ " '^^ — — -, which 

satisfies energy and data constraints, that transmits the same amount of data while having less 
energy expenditure. Hence, we have proved by contradiction that if the rate changes, it must 
increase. 

Proof of Lemma 6 

This proof is similar to the proof of Lemma 5. The main difference is that now the discontinuity 
is in the lower constraint. Then, by following the same procedure we can first show that a rate 
change is suboptimal unless D*^"^\irn) = Dl!^l^{i:l^) and a rate decrease is produced. 
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Proof of Lemma 7 

We know that comes from an energy arrival event, otherwise, the problem would not 
have a feasible solution as stated in Lemma 4. This implies that Dlnax{im) = D^T^^m) and 
= D^^^JC). Then, it is obvious that if Z^^TV^) < D^^^Ji+), an overflow is 
produced. Hence, the proof that = Dinax{im) comes from Lemma 2, where we show 

that the optimal data departure curve satisfies that when an overflow is produced, all the data 
has been transmitted and, hence, the rate must change to zero until the following data arrival. 

Proof of Lemma 8 

This lemma states that two events are produced at time One could look at this lemma as 
the union of Lemmas 5 and 6 at the same time instant, hence, the proof has been already given 
in the aforementioned lemmas. 

Proof of Lemma 9 

As we pointed out, this lemma only applies when the problem has a solution. Consequently, we 
know that Z^f^i(0 = and that Z^^Uc) = Hence, this is the overflow 

problem with the particularity that a data arrival is produced at £m. From Lemma 2, the optimal 
solution minimizes the energy lost due to overflow and, thus, D*''"'\£m) = -D^^x(^m)- However, 
in this situation, nothing can be stated regarding the rate in the following epoch. 

Appendix B 

In this appendix, a technical explanation of the algorithm is given. First, we introduce the 
maximum and minimum rates, a concept required to understand the second and third blocks of 
the algorithm that are presented afterwards. 

A. Maximum and minimum rates 

Let M^^x denote the set that contains the rates obtained by joining the reference point, i.e., 
(0,D(™)(0) = 0)^ with the discontinuities from the left ofD^Ut), i.e., the points {z, D^^^x{z~)), 



'Note that the reference point is always (0, 0) as the origin of coordinates is moved at every iteration. 
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\/z G Zrnax, and such that the obtained curve is feasible for t G (0,2;), where by feasible we 
mean that the curve satisfies all the constraints. Similarly, M^]^ contains the rates obtained from 
joining the reference point with the discontinuities from the right of D^^(t), i.e., the points 
(2, D^^l^{z^)), \Jz G Z^Zin^ ^'^d ^"'^h that the obtained curve is feasible in the same interval. An 
example of this can be seen in Figure 8. 

Let b!^}x denote the infimum of the set 'R^max and R^^i^ refer to the supremum of the set 
M^^j. Let Zmax and z)^^ denote the time instants from which R^}x and -R^^ have been obtained. 
Then, all the rates above Rmax and below -Rf™]i are suboptimal as they would require a rate 

change to transmit the same amount of data.^° Thus, the optimal rate lies within the interval 

roM dM 1 



B. Finish transmission at a constant rate (checkFinish) 

The first step, which is presented in Function 1, checks whether it is possible to transmit all 
bits by using an even power allocation in just one epoch. If it is possible, which implies that 

'"Note that Ul^ax is always greater than -R^^, otherwise, either rI^Jx or -RJ^^ would not be feasible. 
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transmission is ended, the algorithm returns the rate and length of the last epoch, otherwise, it 
returns the strategy or mode that will be used in order to determine the following epoch. 

The function checkFinish first checks whether by transmitting at the maximum feasible rate, 
Rmax, it is possible to transmit all the remaining data -D^] (this is done by the subroutine 
getDatalnCrossing). In case it is not possible, the function returns the mode minEnergy. 
Otherwise, the function finds the time Tq'"^ required to transmit the remaining data -D^] with 
the iteration initial battery Eq"'^^^ Then, it computes the equivalent rate M"^^ and checks whether 
transmitting at this rate is feasible, i.e., the following two conditions are fulfilled: (i.) R^"^'^ < 
Rmax and (ii.) ^('"^ > R^Zir. In case that (i.) is not fulfilled, the function returns the mode minT. 
If (ii.) is not met, it is checked if, by using the following energy arrivals, a feasible curve is 
obtained. Finally, in case both conditions are fulfilled, it is checked whether any energy arrival 
has been produced in the time interval {0,Tq"^^). In case of no arrivals, the algorithm ends and 
the last epoch has been found. In case there is an energy arrival in (0, Tq'"''), the function repeats 
the whole process but now using the initial battery, Eq"^\ and the energy of the first arrival, 
^(m) rpj^.^ process is repeated until (i.) becomes false or a feasible curve is found. 

C. Get rate and length of the next epoch (getEpoch) 

This algorithm's block uses the parameter mode, which is obtained from the function checkFinish 
as presented in Appendix B-B, to compute the rate and length of the following epoch. 

Minimize the total completion time {mode == minT): This strategy is used when both of 
the following conditions are satisfied: (i.) It is possible to finish the transmission at some rate r 
with r < R^max- (ii.) The rate obtained from an even power allocation .R^™^ is not feasible due 
to -R^™^ > R^max- Hence, the objective is to find the rate that allows us to finish transmission as 
soon as possible, without paying attention on saving power, however, without wasting it, either. 
In such case, the rate and duration of the epoch are = R^m}x and irn = zlnJx- 

Minimize the energy expenditure {mode == minEnergy): This strategy is used when, due to 
the constraints, transmission cannot be finished at constant rate, e.g., the rate must be increased 

"Remember that iJg™' = BA{Tm', Tm) — E*{Tm) as summarized in Table I. 
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to satisfy the QoS constraints. Hence, the objective is to save as much energy as possible in order 
to use it when ending the transmission is feasible. Note that in this situation, the problem of 
obtaining the following epoch is similar to the problem presented in [1] and, hence, the solution 
is also similar. The possible data departure curves with constant rate r, i.e., 

= rt, are 

divided in two sets. The first set sIT'' contains all the rates r such that the associated data 
departure curve crosses the constraint DiTax(t) first. Whereas the set contains all the rates 

r such that the associated data departure curve crosses the constraint D^^l^{t) first. Then, the rate 
of the following epoch is determined as the infimum of S^^^^ or, equivalently, the supremum of 
^fllL' ^ i^^Rlax) ^ ^^'^ '•^^ duration of the epoch can be obtained 

as the first time instant such that, rm.im = DlnJx{im) or r^^m = Dl^l{£m)- 

D. Proof of the algorithm optimality 

At each iteration, the algorithm checks whether Lemma 4 is fulfilled. In such a case, there 
is no solution for the problem and the algorithm ends. Otherwise, the algorithm must satisfy 
Lemmas 5-9 at each rate change and Lemma 3 by the end of the transmission. To show that these 
lemmas are satisfied and, hence, the algorithm computes the optimal data departure curve, we 
focus on the three different situations that can occur depending on the constraints of the problem: 
(i.) The algorithm finishes transmission by using an even power allocation among all bits to be 
transmitted, (ii.) minT strategy is used to obtain the epoch, (iii.) The epoch is computed by 
using minEnergy mode. In the following, the optimality of these cases is proved by showing 
that the algorithm-chosen solution satisfies the optimality conditions and that it is unique. 

Part 1 (Even power allocation). When this situation occurs, the algorithm ends transmission 
by transmitting at constant rate. Hence, Lemma 1 is satisfied. No overflow is produced, thus. 
Lemma 2 is satisfied. Lemmas from 5 to 9 do not apply since there are no rate changes. Finally, 
note that Lemma 3 is satisfied, since the last rate tm-i is obtained as the rate that allows the 
transmission of all the data by using all the available energy, moreover, from the properties of 
the function g{-) this rate is unique. 

Part 2 (minT mode). We want to demonstrate that the optimal departure curve transmits at Rlnax 
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during a period of time of zlnJx- Let T'^'") be the total completion time that would be obtained if 
transmitting at rate was feasible, hence, f (™) = D^^l/R^"'\ Similarly, T^nll = D^^l/R^i. 
Note that, since -R^l is a feasible rate, Til^i is an upper bound of the remaining completion 
time, T^™), whereas T^™^ is a lower bound, hence: 

rp{m) ^ rp{m) ^ rp{m) ("yx 
max ■ ^ 

Consider the data departure curve D^^\t) = RtJxt. Note that any other data departure 
curve, D^\t), is suboptimal since, in order to satisfy (7), D^\t) will cross D^;^\t) for some 
ty G (0,Til^l). Hence, at ty, both curves have sent the same amount of data, however, D^{^'\t) 
has consumed less energy. 

Now we must show that, at z'i^ax, the rate increases. Note that by transmitting at R^max instead 
of at some energy has been saved. Then, in the following epoch, the available energy per 

bit is higher and, then, the new rate obtained from an even power allocation among all 

bits in the following epoch fulfills ^('"+^) > > R^ax- Hence, we have proved that at z^}x 
a rate increase is produced. 

Part 3 (minEnergy mode). This mode is used when it is not possible to finish transmission 
at any rate. Note that the algorithm can select three different kind of points, denoted as = 
(irn, D*^"'\im)), for ending the epoch depending on the constraints: 

. vi I = oi^^'Jxii'^} where im is either d'f^ or ef^\ 

. V, I D*'-\i^) = Dit(C) where is either gi™^ or ^ 

. f3 I = DinJxii'J where irn. is ej^^ and an overflow is produced. 

Note that Corollary 1 and Lemma 2 are satisfied for any of the selected points and that Lemma 
3 does not apply since transmission cannot be ended, yet. Hence, we have to prove the following 
three conditions: (i.) If a point such as vi is selected, there will be a rate increase (Lemma 5 or 
Lemma 8 for the case that at two events are produced), (ii.) If a point such as V2 is selected, 
there will be a rate decrease (Lemma 6 or Lemma 8 for the case that at £m two events are 
produced). And (iii.), if a point such as is selected, the rate of the following epoch is zero as 
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far as £3 is not a data arrival event (Lemma 7 and Lemma 9). 

Regarding (i.), the rate of the epoch is = "^"^^ that is the supremum of§i^^^^. Note that, 
at the following iteration, the set ^^^^^ includes all the rates in the interval (0, r^) and the rates 
contained in the interval (r^, + e), for some e > 0. Hence, the rate of the following iteration, 
Tm+i, satisfies that rm+i > r^, therefore, a rate increase is produced. A similar approach can 
be done for (ii.), the rate of the epoch is = that is the infimum of S/jJ^^- At the 

following iteration the set '^^^^^^J contains all the rates in (r^ — e, 00) for some e > 0. Hence, 
the rate of the following iteration, r-m+i, satisfies that r„i+i < r„i and, therefore, there is a rate 
decrease. Finally, in case (iii.), an overflow of the battery is produced. Note that the solution 
chosen by the algorithm satisfies that all the available data at has been transmitted. This 
implies that if non-data arrival is produce at £3 the rate of the following epoch must be zero in 
order to satisfy data causality constraint. 

Therefore, the algorithm computes the optimal solution since it satisfies all the lemmas that 
model the behavior of the optimal solution. 
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Definition 


General notation 


Notation at the m-th algorithm iteration 


Instantaneous power 


P{t) 




Instantaneous rate 


r{t) 




Power-rate function 


<?(•) 




Transmission completion time 


T 


= T - r,„ 


Data departure curve 


D(t) 


L»('")(t) = Z)(t + r^)-D*(r„) 


Optimal data departure curve 


D*{t) 


D*""' (t) = D*{t + r™) - D*{t^) 


Accumulated Data 


DA{t) 


D^2''\t) = DA{t + r^)- D*{r^) 


Minimum Data Departure 


DQos{t) 




Energy expenditure curve 


E{t) 


E^'"^t) = E{t + T,^)-E*{T„,) 


Optimal energy expenditure curve 




E*^"\t) = E*{t + T^)- E*{t^) 


Accumulated Battery 


BA{t;t:,) 


B^^\t) = BA{t + r,„; r™) - E*{t^) 


Minimum Energy Expenditure 


Ejnin (t-, tx ) 




Actual mapping of BA{t;Tm) to data domain 






Effective mapping of BA{t', Tm) to data domain 






Effective mapping of Emin{t',Tm) to data domain 




E'^Eljt) 


Equivalent upper bound on the data domain 






Equivalent lower bound on the data domain 






Discontinuities of the upper bound 




_ u 1 r){'") 1 r)("i) i-t^w 


Discontinuities of the lower bound 






m-th epoch rate 






m-th epoch length 






Beginning m-th epoch 






Data arrival time 


di 


lira) J 

d] ' = di — r,n 


Amount of data in the packet 


A 


d'™' is a relabeling of Di for di > Tm 


Energy arrival time 




(m) 

S-j — Ej Tm 


Amount of energy in the packet 


Ej 


S^""' = B.4(-r„i; Tm) ~ E*{Tm.) 

Ej"^^ is a relabeling of Ej for ej > Tm 


QoS requirement arrival time 


Qk 


(m) 

Ik = Qk - Tm 


Amount of data in the QoS requirement 


Qk 


Q^""' is a relabeling of Qk for > Tm 



Table I 



Summary of the paper notation. 
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Function 1 checkFinish 

Dip^l = D^nJx (oo) > Remaining data 

if {getDataInCrossing{ R^nJxt , D^iTaxit) ) < -D^^t ) then 

return mode = minEnergy, finish = > It is not possible to finish yet. 

else 

for i = : J*^™) — 1 do > J^"^^ is the number of packets with Cj > r„i plus one. 

f/"*) is obtained by solving g~\E/fi'^^) = D^^jf^^ 

Rim) = Dlf^l/f^^"'^ > Even power allocation among all bits 

if R^ZL < R^""^ < R^ax then > D{t) = R^^H is feasible in [0, f}""^] 

§ = {E(-)|j>^,e5.'")G(0,f/'"))} 

if § = then > The algorithm ends and the rate and length of the last epoch are 

returned. 

return r = R^'^l i = f^"'\ finish = 1 
end if 

else if > i^^Tal) then 

return mode = minT, finish = 
end if 
end for 

return mode = minEnergy, finish = 
end if 



